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Combining first-principles density functional theory simulations with IR and Raman experiments, 
we determine the frequency shift of vibrational modes of CO2 when physiadsorbed in the iso- 
structural metal organic framework materials Mg-MOF74 and Zn-MOF74. Surprisingly, we find 
that the resulting change in shift is rather different for these two systems and we elucidate possible 
reasons. We explicitly consider three factors responsible for the frequency shift through physiab- 
sorption, namely (i) the change in the molecule length, (ii) the asymmetric distortion of the CO2 
molecule, and (iii) the direct influence of the metal center. The influence of each factor is evaluated 
separately through different geometry considerations, providing a fundamental understanding of the 
frequency shifts observed experimentally. 
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I. INTRODUCTION 

physiabsorption of small molecules in multi-porous ma- 
terials such as zeolites and metal organic framework 
(MOF) materials has experienced a surge of interest due 
to its potential for hydrogen-storage and gas-separation 
applications^— MOF structures have been widely in- 
vestigated in order to find faster absorption and higher 
storage densities, as well as proper binding energies^— 
In order to design new MOFs with improved properties, 
it is of critical importance to understand the nature of 
the interaction between the absorbed molecule and the 
MOF host. Such understanding can either be gained 
through theory, using first-principles simulations^— or 
through experimental probes, such as infrared (IR) ab- 
sorption and Raman scattering ^ 11 ' 18 ' 20 " — 

Much progress has been made in improving the prop- 
erties of MOFs. For example, it has been shown that 
using unsaturated metal centers, such as MOF74 with 
open Mg, Mn, Zn, Ni, Cu, or HKUST-1 with Cu, 
results in higher absorption density for hydrogen and 
faster absorption at small partial CO2 pressures, the 
latter of which is highly desirable for CO2 capturing 
applications^— ^ ~ 17 ' 19 ~ — It has further been shown that 
iso-structural MOFs with different open metal centers 
can have very different absorption rates at low pressure* 7 - 
In particular, the electronically similar metal centers Mg 
and Zn result in a much faster CO2 uptake rate in Mg- 
MOF74 than in Zn-MOF74 at a pressure smaller than 
0.1 atm* 7 . However, when previous research yields con- 
tradictory results, it becomes difficult to gain further 
insight; e.g. while the frequency shift of the asymmet- 
ric stretch mode of absorbed C0 2 in Mg-MOF74 has 
been reported to be blue shifted in one work using IR 
spectra and B3LYP-D* calculation^ it was reported 



as a red shift in another work using density functional 
theory (DFT) with local density approximation (LDA) 
simulations^ Furthermore, a clear correlation between 
the frequency shifts of the absorbed molecules and other 
absorption properties such as the binding energy or the 
adsorption site is still missing, which makes it difficult to 
directly correlate the observed results with the physical 
nature of the absorption process. 

Van der Waals density functional theory- 3 -^— (vdW- 
DF2) can be used as a very effective tool to under- 
stand the molecule/MOF interactions. Unlike previ- 
ous simulations using LDA or GGA£ where long-range 
van der Waals interactions are not included consis- 
tently, or B3LYP-D*— where the empirical parameters 
are used to incorporate the long-range dispersion terms, 
in our vdW-DF2 method, the exchange-correlation func- 
tional includes the — for these systems so important — 
long-range van der Waals interactions between the MOF 
structure and the physiadsorbed CO2 molecules in a 
seamless^ 2 - and fully self consistent way^ 3 - vdW-DF2 and 
its predecessor vdW-DF have been successfully applied 
to many van der Waals systems^ 4 - ranging from sim- 
ple dimers 3 -^ and physiadsorbed molecules 3 ^ to DNA and 
drug design^ 7 - In particular, it has been demonstrated 
that vdW-DF2 can correctly capture the interaction and 
determine the adsorption sites, binding energy, and vi- 
brational frequencies for small molecules absorbed in dif- 
ferent MOFs£>±i-i£ 

To resolve the contradicting results in literature, and 
to achieve a better understanding of the physics that de- 
termines the frequency shift of the absorbed molecule, 
in this work we combine both theoretical first-principles 
electronic-structure simulation using vdW-DF2 and ex- 
perimental IR and Raman spectroscopy to study the CO2 
absorption in the iso-structural MOFs Mg-MOF74 and 



Zn-MOF74. In the first step, we use our experimental 
probes to determine the absorption behavior and the vi- 
brational frequencies of CO2 absorbed in these MOFs. 
Then, in the next step, we use vdW-DF2 to calculate 
the corresponding shifts and determine the mechanisms 
that cause them. Unlike in experiments, first-principles 
simulations allow us to artificially freeze different de- 
grees of freedom of the system, enabling us to understand 
the importance of different physical contributions to the 
CO2 frequency shifts. By analyzing the CO2 asymmetric 
stretch frequency under different geometries, we identify 
three different factors determining the frequency shift of 
physiadsorbed CO2, namely, the length of the molecule, 
the asymmetric distortion, and the metal center. 



II. METHOD 

A. Metal organic framework synthesis 

Mg-MOF74^ Zn-MOF74,^ Co-MOF74,-22 and Ni- 
MOF74£ were synthesized according to procedures de- 
scribed in the literature^^ 

Mg-MOF74: 2,5-Dihydroxyterephthalic acid (H 2 DHB 
DC) (99 mg, 0.5 mmol) and Mg(N0 3 )2-6H 2 
(257 mg, 1.0 mmol) were dissolved in the mixture of 
THF (7 mL), 1M NaOH solution (2 mL), and water 
(3 mL) with stirring. The mixture was then sealed 
in a Teflon-lined autoclave and heated at 110°C for 
3 days. The product was collected by filtration as 
a light-yellow substance. Yield: 115 mg, 83%. 

Zn-MOF74: The preparation of Zn-MOF74 was similar 
to that of Mg-MOF74 except that Zn(N0 3 ) 2 -6H 2 
(298 mg, 1.0 mmol) was used instead of Mg(N0 3 )2- 
6H 2 0. Yield: 160 mg, 87%. 

Co-MOF74: H 2 DHBDC (150 mg, 0.75 mmol) and 
[Co(N0 3 ) 2 ]-6H 2 (186 mg, 0.75 mmol) were dis- 
solved in 15 mL of THF-H 2 solution (50:50, v:v) 
with stirring. The mixture was transferred to a 
Teflon-lined autoclave, which was then sealed and 
heated at 110°C for 3 days. Brown-red rod-shape 
crystals were isolated by filtration and dried under 
vacuum. Yield: 130 mg, 50%. 

Ni-MOF74: A mixture of H 2 DHBDC (60 mg, 
0.3 mmol), [Ni(N0 3 ) 2 ]-6H 2 (174 mg, 0.6 mmol), 
DMF (9 mL), and H 2 (1 mL) was transferred to 
a Teflon-lined autoclave and heated in an oven at 
100°C for 3 days. Brown crystalline powder was 
collected by filtration and dried under vacuum. 
Yield: 75 mg, 72%. 

All as-synthesized materials were exchanged with fresh 
methanol four times in a duration of 4 days, followed 
by drying in an vacuum oven at room temperature, and 
annealing at 480 K overnight under high vacuum before 
spectroscopic measurements. 



B. IR and Raman spectroscopy 

IR absorption spectroscopy of CO2 absorption in the 
MOFs was performed in transmission at room temper- 
ature using a liquid-N 2 -cooled InSb detector. Approxi- 
mately 12 mg of MOF powder was pressed on a KBr sup- 
port and mounted in a high-temperature high-pressure 
cell (Specac product P/N 5850c) and heated to 200°C in 
vacuum (100 mTorr) overnight for complete dcsolvation. 
MOF74 samples were activated by solvent exchange in 
methanol and drying in vacuum at room temperature. 
Subtraction of the gas phase C0 2 contribution to the IR 
spectra was performed as described in Ref. [ill- 
Raman spectroscopy measurements were performed 
using a solid state 532 nm laser. The activated sam- 
ple was loaded into a Linkam FTIR cooling/heating 
stage, and the sample was heated to 120°C in vacuum 
(900 mTorr) for complete dehydration. A laser power 
of 0.113-1.23 mWatt was used to avoid sample burning 
from the laser. 



C. First-principles calculations 

For our first-principles calculations we used DFT as 
implemented in ABINIT^S. utilizing vdW-DF2 to de- 
scribe exchange and correlation effects^ Troulier-Martin 
type norm-conserving pscudopotcntials and a plane- wave 
basis are used,— where the Zn 3d semicore electrons are 
also considered as valence electrons. To ensure a full 
convergence of the structure and energy, a kinetic en- 
ergy cutoff of 45 Hartree is used for the plane- wave ba- 
sis. For structural relaxation, we start from the experi- 
mental atomic positions and relax them using vdW-DF2 
until the force on each atom is less than 0.05 eV/A. The 
unit cell parameters are fixed to the experimental val- 
ues, where for the hexagonal unit cell a = 25.881 Aand 
c = 6.8789 A for Mg-MOF74^ and a = 25.887 Aand 
c = 6.816 A for Zn-MOF74£ respectively. 

Due to the complex interaction between CO2 and 
MOF74, many local minimum energy sites exist, mak- 
ing it difficult to identify the low-energy adsorption site, 
which corresponds to the low pressure absorption. In- 
deed, multiple adsorption sites and binding energies are 
determined from our simulations. In this work, we per- 
formed calculations with multiple initial positions and 
orientations of CO2 relative to the MOF74, and choose 
the site with the lowest total energy as the one for our 
analysis. To guarantee that CO2 is absorbed at the 
equivalent lowest energy site in both Mg-MOF74 and Zn- 
MOF74, we start the atomic position relaxation of CO2 
in Mg-MOF74 with the coordinates of the relaxed CO2 
in Zn-MOF74, which is possible due to the isostructure 
and symmetry of the two MOF systems. This ensures 
that we are comparing the effect of the open metal site 
on CO2 in the two different systems on the same footing. 

To calculate the frequency shifts we use a frozen- 
phonon approach in the CO2 molecules, where the 




FIG. 1. (Color online)Illustration figure of C0 2 absorbed in Zn-MOF74. 
Zn-MOF74 and Mg-MOF74 can be found in Table [Q 



Calculated values of di, di, d3, and 9 for CO2 in 



TABLE I. Calculated binding energy AE (kj/mol), angle 
8 (deg), and various distances di and I (A) (where I — di+dz) 
of the CO2 molecule physiadsorbed in the MOF ; see Fig. Q] 
for further details. The calculated free CO2 parameters are 
d = 1.1630 A and l Q = 2.3260 A. 



MOF 


AE 
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di 


di — do 


di — do 


l-lo 


Mg-MOF74 
Zn-MOF74 


35.4 
26.9 


120.8 
116.4 


2.53 
2.69 


+0.0062 
+0.0048 


-0.0059 
-0.0039 


+0.0003 
+0.0009 



MOF74 atoms are kept fixed, while each of the atoms 
in the CO2 is distorted in ±dx, ±dy, ±dz directions by 
small distortions Ar — 0.02 A, to calculate the force 
on each atom. The symmetrized force matrix 2F(Ar) = 
F(Ar) — F(—Ar) is constructed for only the atoms within 
the CO2 molecule, based on the approximation that the 
interaction between the CO2 molecule and the MOF is 
weak and thus the effect of the vibration of the MOF74 
atoms on the CO2 frequencies are negligible. However, 
notice that, while the MOF74 atoms arc kept fixed, the 
van der Waals forces experienced by the CO2 due to the 
presence of the MOF74 are included by our vdW-DF2 
calculation. To evaluate the effects of the surrounding 
MOF lattice, tests have been performed to allow vibra- 
tions of several MOF74 atoms in the vicinity of the ab- 
sorbed CO2 molecule; frozen phonon calculations with 
such an extended force matrix show no observable effect 
on the CO2 frequencies. In general, our studies show that 
calculated frequencies are converged to within less than 
1 cm -1 . 



III. RESULTS AND DISCUSSION 
A. Structure and binding energy 

The structure of CO2 absorbed in MOF74 is illus- 
trated in Fig. [1] In this work, we consider only CO2 ab- 
sorption under low pressure smaller than 7 Torr, where 
the open metal site in MOF74 is far from fully occu- 
pied. In our first-principles simulations one CO2 per six 
metal sites absorption is considered, corresponding to 



the low-loading situation observed experimentally. For 
the binding energy of CO2 absorbed in Mg-MOF74, we 
find 35.4 kJ/mol, while it binds somewhat weaker in Zn- 
MOF74 with 26.9 kJ/mol. Comparing the distance of 
the CO2 molecules with the open metal site, we find a 
metal-oxygen distance (d\ in Fig.[lJ of 2.53 A for Mg ver- 
sus 2.69 A for Zn. In other words, CO2 binds stronger 
and closer to the metal center in Mg-MOF74 than it does 
in Zn-MOF74. At the same time, similar metal-C02 an- 
gles ensure that the CO2 is absorbed at equivalent sites 
in both MOFs, with slightly different metal-C02 angles 
of 120.8° and 116.4° in Mg-MOF74 and Zn-MOF74, re- 
spectively. 

Further analyzing the structure of the absorbed CO2 
molecule, we see that the C=0 bond closer to the metal 
site (^2 in Fig. [1]) is elongated in both systems with a 
value of 1.1692 A in Mg-MOF74 and 1.1678 A in Zn- 
MOF74; in free C0 2 it is 1.1630 A (denoted as d in Ta- 
ble HJ. On the other hand, the C=0 bond farther away 
from the metal center (g^ in Fig. [IJ are both shortened, 
with a value 1.1571 A in Mg-MOF74 and 1.1591 A in Zn- 
MOF74. A summary of distortions is shown in Table [J] 
This distortion can be understood intuitively via the in- 
teraction between the CO2 and the metal center, where 
the attraction from the metal center weakens (and thus 
elongates) the nearby C=0 bond, shortening the remain- 
ing C=0 bond. By comparing with the free CO2 value of 
1.1630 A, we see that the CO2 asymmetric distortion in 
Mg-MOF74 is stronger than in Zn-MOF74. Summing up 
the two C=0 bonds yields the overall length of the CO2 
molecule, which shows an elongation of +0.0003 A and 
+0.0009 A for Mg-MOF74 and Zn-MOF74 absorption, 
as shown in Table HI It has been reported previously that 
CO2 molecules absorbed in the MOF structure might ex- 
perience some nonlinear distortion, i.e., an O-C-0 angle 
differing from 180°—. In this work, we also alow nonlinear 
distortion of the CO2 molecule, and we find that the re- 
laxed CO2 molecule are only slightly bent after adsorbed 
within MOF74, featuring an O-C-0 angle of 179.25° in 
Zn-MOF74 and 178.97° in Mg-MOF74. 

As mentioned above, multiple local minimum adsorp- 
tion sites arc obtained during the relaxation of CO2 
within the MOF74, leading to different mctal-C02 dis- 



tances and binding energies. For instance, in Zn-MOF74 
positions with similar location but Zn-C02 distances of 
2.96 A and 3.09 A arc found with binding energies of 
25.6 kJ/mol and 24.8 kJ/mol, respectively. In addition 
to these sites similar to the global minimum metal site, 
we also find a secondary non-metal site for CO2 absorp- 
tion, which lies between two neighboring equivalent metal 
sites along the direction of the MOF74 pore, with a dis- 
tance of 3.99 A to the nearest metal Zn site and a binding 
energy of 21.8 kJ/mol — clearly much smaller than at the 
metal site. The distance between two equivalent metal 
sites along the pore direction in MOF74 is the c length 
of the hexagonal unit cell, which is approximately 6.8 A, 
as mentioned above. It is thus possible that at high CO2 
pressure, when all the metal sites are occupied (1 CO2 
per 1 metal), more CO2 can still be absorbed within the 
MOF74 occupying the secondary site, with a CO2-CO2 
distance of approximately 3.4 A along the pore direction, 
thus resulting in a higher (2 CO2 per 1 metal) storage 
density of C0 2 within MOF74. 



B. CO2 frequencies from experiment 

IR absorption spectra of CO2 absorbed in Mg-MOF74 
and Zn-MOF74 are shown in Fig. [21 depicting a main IR 
absorption band attributed to the asymmetric stretch of 
C0 2 at 2338 and 2352 cm" 1 for Zn-MOF74 and Mg- 
MOF74. The integrated areas of the IR band of the 
asymmetric CO2 stretch mode can serve as an indica- 
tion of the relative amount of CO2 absorbed in MOF74. 
Figure [3] summarizes the integrated area as a function 
of C0 2 pressure for both Mg-MOF74 and Zn-MOF74, 
showing that under the same pressure, the amount of 
absorbed CO2 is more for the case of Mg-MOF74 than 
it is for Zn-MOF74. These results are consistent with 
isotherm measurements of CO2 reported in Ref. |7J. Sim- 
ilar IR measurements are also performed on Co-MOF74 
and Ni-MOF74, showing shifts of the asymmetric CO2 
stretch (-2340 cm -1 ) similar to that in Zn-MOF74. 

In Fig. [2 the shoulder peaks at 2325 and 2341 cm" 1 arc 
attributed to the combination mode of the stretch mode 
and the two non-degenerate bending modes (denoted as 
Si and 8 2 in Fig. [2]). Experimentally, it is difficult to 
identify the bending modes in the low frequency range 
because of their weak intensity as compared to the MOF 
vibrations. Therefore, we could only resolve one of these 
non-degenerate bending modes at —658 cm -1 . The sim- 
ilarity of the bending mode for the Mg and Zn cases is 
consistent with the close calculated values for both sys- 
tems shown in Tabic [TTJ 

Raman spectroscopy measurements were preformed to 
find the experimental value of the symmetric stretch 
mode of absorbed C0 2 for the Co-MOF74 and Ni- 
MOF74. The symmetric stretch was found to be ap- 
proximately -1380 and -1382 cm -1 for Co-MOF74 and 
N1-MOF74, respectively. Because of the strong MOF Ra- 
man modes, the symmetric stretch for Mg-MOF74 and 
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FIG. 2. IR absorption spectra of CO2 absorbed into Zn- 
MOF74 (top) and in Mg-MOF74 (bottom) at changing CO2 
pressure (1-6 Torr). 



Zn-MOF74 is not available. A summary of the IR and 
Raman data for the frequencies of the absorbed CO2 
within MOF74 can be found in Table UU 



C. Analyzing the frequency shift 

When wc consider the frequency shifts from the free 
CO2 situation, comparing with 2349 cm -1 for free CO2, 
the asymmetric stretch for CO2 absorbed in MOF74 
has been shifted by +3 and -11 cm -1 for Mg and Zn- 
MOF74, respectively Our first-principles frozen-phonon 
calculations are in good agreement with this result and 
we find the CO2 asymmetric stretch mode experiences a 
slight redshift of approximately -0.5 cm" 1 in Mg-MOF74 
and -8.1 cm" 1 in Zn-MOF74. Furthermore, the bend- 
ing modes calculated also experience redshifts of about 
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Zn-MOF74 as a function of CO2 pressure. 



-10 cm- 1 in both Mg-MOF74 and Zn-MOF74, which is 
in excellent agreement with the experimentally observed 
-9 cm -1 rcdshifts in these systems. The calculated re- 
sults for all the CO2 modes are summarized in Tabic |TTJ 

We have shown above that CO2 is binding stronger 
in Mg-MOF74, with a closer mctal-C02 distance and 
larger molecule distortions. It is thus puzzling to ob- 
serve that the ^3 frequency shift for CO2 in Mg-MOF74 is 
much smaller compared with that in Zn-MOF74. There 
have been many experiments and theoretical simulations 
studying the frequency shifts of small molecules absorbed 
in MOFs, such as H 2 , C0 2 , CO, etc^l^rM. It is well 
known that the frequency shift of small molecule vibra- 
tions has no obvious correlations with either the bind- 
ing energy or the adsorption site of the molecule within 
the MOF. In our example here with iso-structural Mg- 
MOF74 and Zn-MOF74, this behavior is even more obvi- 
ous, as Mg and Zn have similar electronic structures: the 
metal centers have valence electrons of 3s and 4s, respec- 
tively, except that Zn has the additional fully occupied 
3d orbital electrons. Thus the question arises as to what 
the physical reason is that determines the frequency shift 
of the small molecules absorbed in the MOF. 

Several reasons may contribute to the frequency shift. 
For instance, by analyzing the geometries of the absorbed 
CO2, we noticed that the absorbed CO2 molecules have 
been distorted from their free molecule geometry, exhibit- 
ing an off-center shift of the carbon atom, as well as an 
elongated overall length of the molecule, as summarized 
in Table [T^The slight nonlinear distortion of the CO2 
molecule is ignored here in this analysis). The change 
in the length of the molecule and the asymmetric distor- 
tion can both affect the vibrational frequencies. Beyond 
the effect of the change in the molecule geometry, the 
nearby open metal center might also play a direct role in 
the frequency shift, by attracting or repelling the nearby 
oxygen atom in the CO2 molecules during the vibration. 



TABLE II. Vibrational frequencies (cm- 1 ) of CO2 physiad- 
sorbed in MOF74. Higher accuracy is kept for the calculated 
vz for further detailed analysis. 



system 



sym. v\ 



bend vi 



asym. 1/3 



exp. 



free CG^ 


1388 


667 


2349 


Mg-MOF74 


— 


658 


2352 


Zn-MOF74 


— 


658 


2338 


Co-MOF74 


1380 


— 


2340 


Ni-MOF74 


1382 


— 


2340 





free C0 2 


1298 


646, 639 


2288.5 


calc. 


Mg-MOF74 


1294 


636, 630 


2288.0 




Zn-MOF74 


1296 


637, 632 


2280.4 



Taken from Ref. Wl 



To understand these different contributions, in the fol- 
lowing we analyze the influences of each of these factors 
on the frequency of the CO2 vibrations separately and 
summarize their contributions in Tabic IFlTI 



1. Change in molecule length 

We first analyze the effect of the changing molecule 
length. To do this, wc perform frozen-phonon calcula- 
tions on the free CO2 molecules fixed to the lengths of 
the CO2 absorbed in Mg-MOF (shown in the last col- 
umn of Table H|, while keeping the carbon atom at the 
center of the molecule. With this geometry, the CO2 
asymmetric stretch is shifted by —1.6 cm . A simi- 
lar calculation by setting the CO2 molecule length to 
the value of that absorbed in Zn-MOF74 yields a red- 
shift in v-i of -3.7 cm -1 . Note that overall the CO2 
in Mg-MOF74 is elongated by 0.0003 A, while that in 
Zn-MOF74 is elongated by 0.0009 A. The C0 2 with the 
longer length (in Zn-MOF74) has a larger redshift of ap- 
proximately -2.1 cm" 1 than that in Mg-MOF74. That 
is, a longer molecule has more redshift, as suggested by 
common sense. 



2. Asymmetric distortion of the molecule 

Next, we take the asymmetric effect into considera- 
tion by placing the CO2 at exactly the same geometries 
as they are in the two MOFs, while removing the sur- 
rounding MOFs. The frequency shift thus comes from 
the change in both the length of the CO2 molecule and 
the asymmetric shift of the carbon atom. Our frozen- 
phonon calculations give an overall frequency shift of -0.5 
and -3.0 cm -1 for the 1/3 mode of C0 2 from Mg-MOF74 
and Zn-MOF74, respectively. Subtracting the previous 
results of only considering the length of the molecule, 
the asymmetric distortion of the carbon atom causes a 
slight blue shift in the asymmetric stretch with a value 
of 1.1 and 0.7 cm- 1 in Mg-MOF74 and Zn-MOF74. 



To confirm these results, we also perform calculations 
using another setup, where we keep the optimized CO2 
within the MOF and shift the carbon atom to the center 
of the CO2 molecule. In this way, we include the effect 
of the MOF environment and the length effect. The only 
difference between this new setup and the optimized CO2 
in MOFs is that the carbon atom off-center asymmetric 
effects are eliminated. Our simulations show that the 
CO2 in such a geometry has a frequency shift of —1.7 and 
-8.7 cm' 1 in Mg-MOF74 and Zn-MOF74. Comparing 
these results with the frequency shifts for relaxed CO2 in 
the same MOFs, we find similar differences of approxi- 
mately 1.2 and 0.6 cm' 1 for the case of Mg-MOF74 and 
Zn-MOF74, which are consistent with our previous re- 
sults. It is thus clear that the carbon off-center distortion 
indeed causes slight blue shifts in the 1/3 frequency. This 
result can also be understood in the following way: by 
shifting the carbon atom off-center, the CO2 atom is dis- 
torted asymmetrically in a similar pattern as the asym- 
metric stretch mode, which thus favors the asymmetric 
stretch. As shown in Table HI the asymmetric distortion 
in Mg-MOF74 is stronger than that in Zn-MOF74, which 
is consistent with a larger blueshift effect in the CO2 ^3 
frequency 



3. Effect of the metal center 

By comparing with the results for optimized CO2 in 
the MOFs, the previous results of the effects of the CO2 
molecule geometry change on the v% frequency shifts yield 
the influence of the metal center on the V3 frequencies. A 
direct comparison gives a frequency shift of -5.1 cm' 1 for 
the Zn-MOF74 and cm" 1 for Mg-MOF74. This result 
is quite a surprise, claiming that the open Zn center has 
a strong redshift influence on the ^3 frequency while the 
Mg center has no effect at all. 

To confirm this result, we place the undistorted CO2 
molecule at the same position and angle of the adsorp- 
tion site, with the position of the oxygen atom near the 
metal center fixed, while the positions of the carbon and 
the other oxygen atom shifted slightly along the line to 
achieve the free CO2 geometry. The frequency shift of 
the CO2 asymmetric stretch thus results mostly from 
the direct effect of the metal center. With this geom- 
etry, our calculation shows that, while the asymmetric 
stretch of CO2 in Zn-MOF74 is shifted by approximately 
-5.0 cm -1 , the one in Mg-MOF74 has a negligible shift 
of about -0.6 cm" 1 , confirming our previous results. In 
other words, the fact that the oxygen atom in CO2 is 
being close to the open metal center has a significantly 
different effect depending on the metal atoms. While 
the Zn atom affects the frequency strongly, the Mg atom 
has almost no effect at all. This result is quite striking, 
since the metal centers Mg and Zn have a very similar 
valence electronic structure with 3s and 4s electrons as 
the outermost valence states. The result thus shows that 
fully occupied scmicorc 3d electrons in Zn have an im- 



TABLE III. Frequency shifts of V3 (cm 1 ) for different ge- 
ometries. See main text for details about these geometry. 



effect 


Mg-MOF74 


Zn-MOF74 


length 


-1.6 


-3.7 


length+asym 


-0.5 


-3.0 


length+metal 


-1.7 


-8.7 


metal 


-0.6 


-5.0 


overall 


-0.5 


-8.1 




FIG. 4. (Color online)Charge density differences Ap — 
Pmof74,+co 2 — PMOF74 — pco 2 from before and after CO2 
adsorbed into Zn-MOF74 (top) and in Mg-MOF74 (bottom). 



portant effect in the interaction with the absorbed CO2 
molecules. In fact, in Co-MOF74 and M-MOF74 where 
3d electrons are also present in the metal center, experi- 
mentally we observed that the asymmetric stretch of CO2 
has been similarly red-shifted by -10 cm" 1 as that in Zn- 
MOF74, as summarized in Table [TT1 . indicating that the 
3d orbitals are indeed playing a similar role for the CO2 
frequencies in these MOF74 systems. 

The contributions of these three effects to the 1/3 fre- 
quency shifts are summarized in Table [TTT1 



D. Microscopic insight and implication of bonding 

In order to gain more insight understanding for the dif- 
ferences in the interaction between CO2 and the two dif- 
ferent MOF74s, we investigate the charge density differ- 
ences before and after CO2 adsorbed inside the MOF74. 
To do this, we calculate the charge densities of the 
CO2 and MOF74 separately, by removing CO2 from the 
MOF74, while keeping the atomic positions. The charge 



density difference is defined by 

A/? = PMOF74+C02 ~ PMOF74 ~ PCOi ■ 

The charge density difference Ap thus illustrates the ef- 
fect of the interaction by placing the CO2 inside the 
MOF74. The obtained charge density differences are 
plotted in Fig. |U 

Not surprisingly, by introducing CO2 inside the 
MOF74, the electrons in CO2 are attracted toward the 
metal side, featuring electron deficiency in the far end 
(light blue color) and electron gaining near the metal 
site (yellow color) . This can be understood via the attrac- 
tion from the positively charged metal centers, similar for 
both Mg and Zn centers. However, one can also observe 
differences between the Zn-MOF74 and Mg-MOF74 sys- 
tems, where the electron transfers more in the later than 
the former. As we look at the metal center, we can ob- 
serve that near Zn atom, the electrons shift slightly away 
from the CO2 adsorption position. This charge density 
change near the Zn centers are missing from those Mg 
center, clearly showing that it is an effect of the d-orbitals 
of the Zn atoms. It is thus clear that the presence of 
the (i-orbitals prevents the charge density transfer within 
the adsorbed CO2 molecules, leading to smaller charge 
transfer as compared with those near Mg center. As a 
result, the bonding energy between CO2 and Zn-MOF74 
is smaller than that of C0 2 in Mg-MOF74. On the other 
hand, this clear effect of the Zn <i-orbital effectively af- 
fect the IR frequency shift of the adsorbed CO2 molecule, 
leading to the differences between the two systems as we 
discussed above. 



E. Discussion 

From the previous analysis, it is now clear that the z/3 
frequency shifts of the CO2 absorbed in MOF74 can be 
understood by the three distinct contributions, namely, 
the change of the molecule length, the off-center asym- 
metric distortion, and the direct effect of the open metal 
site. The interaction between CO2 and the MOF causes 
distortions in the molecule geometry, which affect the 
vibration modes. The absorbed CO2 molecules experi- 
ence a stronger asymmetric distortion in Mg-MOF74, as 
shown in Table U which is consistent with the larger 
bonding energy calculated, as well as the larger inte- 
grated area for the asymmetric stretch IR peak. However, 
such an asymmetric distortion only has a small effect on 
the V3 frequency shifts. On the other hand, although 
CO2 is less asymmetrically distorted in Zn-MOF74, it ex- 
periences a larger elongation in the overall length, which 
affects the z/3 frequency shift more strongly. In addition, 
the nearby open metal site can play a quite different role 
in affecting the CO2 vibrations, where in this work, the 
Zn affects strongly while Mg has barely any effect on the 
1/3 frequency shift. The differences between these two 
MOF74 can be intuitively understood by the presence 
of the (i-orbitals within the Zn atoms and the missing 



of those within Mg centers, as illustrated in the charge 
density difference plots. 

The results of this work provide insight to the factors 
that determine the frequency shifts of the absorbed CO2 
in MOF, helping to understand the puzzling frequency 
shifts observed experimentally. More importantly, the 
analysis method of this work can serve as a new way to 
understand the more widely examined molecule-MOF in- 
teractions and frequency shifts. However, one must keep 
in mind that frequency shifts obtained through such fixed 
geometries and environments reflect the influence of dif- 
ferent factors on the force matrix and can only give an 
estimation of the influence of certain factors. In reality, 
originating in the molecule-MOF interaction, all three 
factors are closely connected intrinsically and it is im- 
possible to exactly separate these different effects. 



IV. SUMMARY 

In this work, we analyzed the physics determining the 
asymmetric frequency shift of the CO2 molecules phy si- 
adsorbed in MOFs. Our specific findings are summarized 
as follows: (i)first-principles vdW-DF2 simulations de- 
termine that the C02's have a closer distance to the Mg 
center and a larger binding energy within Mg-MOF74 
comparing with those in Zn-MOF74. (ii) Contrary to 
our intuition, and despite the isostructure and the sim- 
ilarity of the open metals Mg and Zn, the asymmetric 
stretch frequency of physiadsorbed CO2 has been shifted 
stronger in Zn-MOF74 (-11 cm" 1 by IR and -8.1 cm" 1 
by theory) than that in Mg-MOF74 (+3 cm" 1 by IR 
and -0.5 cm" 1 by theory) .(iii) By comparing the re- 
sponse in two isostructure MOFs, namely Zn-MOF74 and 
Mg-MOF74, we identified the three most important fac- 
tors contributing to the frequency shifts: the elongated 
CO2 molecule, the off-center asymmetric distortion of 
the carbon atoms, and the effect of the metal center, 
(iv) The asymmetric stretch frequency is very sensitive 
to the overall length of the CO2 molecule. Absorbed in 
the MOF, the CO2 molecules are elongated, which leads 
to a redshift in the frequency. This elongation effect and 
resulting redshift are more significant for CO2 absorbed 
in Zn-MOF74 compared with those in Mg-MOF74. (v) 
The slight off-center asymmetric distortion, on the other 
hand, favors the asymmetric stretch and causes a slight 
blueshift in the frequency, (vi) Aside from changing 
the geometries of the CO2 molecule (i.e. elongating the 
molecule, causing off-center asymmetric distortion of the 
carbon atom) and depending on the species of the open 
metal site, the direct interaction of the oxygen atom with 
the metal center can have very different effects on the fre- 
quency of the asymmetric stretch, where the Zn center 
leads to a redshift of about -5 cm -1 and the Mg center 
has a negligible effect on the frequency, (vii) The ob- 
served different effects of the Zn-MOF74 and Mg-MOF74 
can be understood by the presence of <i-orbital electrons 
in the Zn-MOF74. 
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